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ABSTRACT 



We report on a timing of the eclipse arrival times of the low mass X-ray binary and X-ray pulsar 2A 1822-371 performed 
using all available observations of the Proportional Counter Array on board the Rossi X-ray Timing Explorer, XMM- 
Newton pn, and Chandra. These observations span the years from 1996 to 2008. Combining these eclipse arrival time 
measurements with those already available covering the period from 1977 to 1996, we obtain an orbital solution valid for 
more than thirty years. The time delays calculated with respect to a constant orbital period model show a clear parabolic 
trend, implying that the orbital period in this source constantly increases with time at a rate Porb = 1.50(7) x 10"^'^ 
s/s. This is 3 orders of magnitude larger than what is expected from conservative mass transfer driven by magnetic 
braking and gravitational radiation. From the conservation of the angular momentum of the system we find that to 
explain the high and positive value of the orbital period derivative the mass transfer rate must not be less than 3 times 
the Eddington limit for a neutron star, suggesting that the mass transfer has to be partially non-conservative. Under 
the hypothesis that the neutron star accretes at the Eddington limit we find a consistent solution in which at least 70% 
of the transferred mass has to be expelled from the system. 
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1. Introduction 

The source 2A 1822-371 is a well-known low mass X- 
ray binary (hereafter LMXB) se en almost edge-on wit h 
an inclination angle of i ^ 85° (|Hellier fc MasonI Il989l) . 
The observed average unabsorbed flux of the source in the 
0.1 — 100 keV energ: y range is ~ 1.5 x 10~^ ergs cm~^ 



s~" (jlaria et al.ll2001[ ). This corresponds to an unabsorbed 
luminosity of ~ 1.2 x 10^^ ergs s~ ^, adopting a distance 
of 2.5 k pc (iMason fc Cordovalll982l ). However, it has been 
noted (jParmar et al. 20001) that the mean ratio of the X- 
ray over optical luminosity, Lx/Lopt, for 2 A 1822-371 is 
about 20, while t he average value for LMXBs is about 500 
()van Paradiis fc M cClintock 1994). This would imply an 
' unobscured X-ray luminosity as high as 3 x 10'^'' ergs/s for 
, the assumed distance of 2.5 kpc. The apparent low lumi- 
nosity of the source has therefore to be ascribed to the high 
inclination of the system with respect to the line of sight. 
Indeed, the light curve of 2A 1822-371 shows both dips and 
eclipses of the X-ray source by the companion star. The 
partial nature of the eclipse indicates that the X-ray emit- 
ting region is extended and that the observe d X-rays are 
scattered in an accretion disk corona (ADC, iWhite et al.l 
(|1981h ). The X -ray light curve shows clear signs of orbital 
modulation with a binary orbital period of 5.57 h. This X- 
ray modulation is probably caused by the obscuration of 
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the ADC by the thick rim of an accretion disk. The or- 
bital period has been measured from eclipse timing to in- 
crease gradually (jHellier et al.lll990[ ). IParmar et alJ feOOO") 
gave the best ephemeris of this source before this work. In 
particular they found a significant positive orbital period 
derivative of Porb = 1-78 x 10"^° s/s. 

iJonker Sz van der KlisI (|200lh reported on the discovery 
of 0.59 s X-ray pulsations in this source in an RXTE ob- 
servation performed in 1998. The timing analysis of the 
pulse arrival times indicates a circular orbit with an ec- 
centricity e < 0.03 (95% c.l.) and an asinz for the neu- 
tron star of 1.006(5) It-s, implying a mass function of 
(2.03 ± 0.03) X 10^2 Mq. The comparison between the 
pulse period measured by RXTE in 1996 and 1998 also 
indicates that the neutron star in this system is spin- 
ning up at a rate of P = (-2.85 ± 0.04) x lO^^^ 
iJonker fc van der KlisI (|200lD inferred a bolometric X-ray 
luminosity of about (2 — 4) x 10'^'' ergs/s assuming a mag- 
netic field of (1 — 5) X 10^^ Gauss. From spectroscopic mea- 
surements of the r adial velocity curve of the companion, 
IJonker et al.l (|2003l ) derived a lower limit to the mass of the 
neutron star and to that of the companion star of 0.97±0.24 
and 0.33 ±0.05 Mq, respectively (1 a, including uncertain- 
ties in the inclination), and an accurate estimate of the 
system inclination angle, i = 82°. 5. 

In this paper we report on the analysis of X-ray ob- 
servations of 2A 1822-371 performed from 1996 to 2008 
by RXTE, XMM- Newton, and Chandra with the aim to 
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derive eclipse arrival times and to improve the orbital 
ephemeris. We confirm with higher precision and over a 
m uch larger time s pan ( about 31 years) the ephemeris found 
bv iParmar et al.l (j2000r) . In particular we find that the or- 
bital period derivative has remained constant during the 
last 30 years. Finally we discuss the implications of a high 
and positive value of the orbital period derivative on the 
mass transfer rate and secular evolution of this source. 

2. Timing analysis and results 

We analysed all available X-ray observations of 2A 1822- 
371 performed over the period from 1996 to 2008. In partic- 
ular we used observations from the PCA on board RXTE 
performed in 1996 (P10115), 1998 (P30060), 2001 (P50048), 
2002 (P70036), 2002-2003 (P70037), one observation from 
XMM-Newton performed in 2001 (Obs ID: 0111230101 and 
0111230201), and two Chandra observations performed in 
2000 (Obs ID: 671) and in 2008 (Obs ID: 9076 and 9858), 
respectively. The arrival times of all events were referred to 
the solar system barycenter, using as the best estimate for 
the source coordinates those derived from the 2008 Chandra 
observations (RA: 18 25 46.81, DEC: -37 06 18.5, uncer- 
tainty: 0.6"). 

The typical eclipse duration is around 2.2 ks, which cor- 
responds to 10% of the binary orbital period. In order to 
improve the statistics for the measure of the eclipse epochs 
and to have the possibility of fitting a complete orbital light 
curve we decided to perform a folding of these data using 
the known binary orbital period of the source, after verify- 
ing that this folding does not affect the results reported here 
in any case. Folding the data is not an important issue for 
the two Chandra observations and the XMM observation, 
where just one or two consecutive eclipses are observed. But 
it is important for the RXTE observations, because these 
are short and sparse, and also because the RXTE observa- 
tions are continuously interrupted by the Earth occultation 
at every RXTE orbit (lasting approximately 1.5 h). In this 
case the folding is required to sample a complete orbital 
light curve from the source, because this is important for a 
meaningful fitting of the eclipse. For each of these observa- 
tions we hence folded the data using the local orbital period 
as der ived from the ephemeris published by iParmar et al.l 
(|2000ll . The 2002-2003 RXTE dataset (P70036 and P70037) 
was long enough and we decided to divide it into the follow- 
ing four periods: i) 2002 June 7-10, ii) 2002 August 2-18, iii) 
2002 September 2-30, and iv) 2003 August 31 - September 

3. In this way we obtained a total of 10 orbital light curves 
in which the eclipses were clearly visible (see Table [1] for 
details on the used observations). 

We then fitted these orbital light curves to derive eclipse 
arrival times with the procedure described below. Because 
the eclipses are asymmetrical and partial, the exact eclipse 
centroid times crucially depend on the model adopted to 
describe their shape as well as the variable continuum they 
are superimposed on. In order to be conservative in our 
estimates, we then decided to fit the folded light curves 
using 10 diffe r ent ni odels. The first model is that used by 
IParmar et al.l (|2000D consisting of a Gaussian and a con- 
stant fitted on a phase interval of 0.1 around the eclipse. 
The second and third models consist again of a Gaussian 
and a constant plus a linear term (second model) and a 
linear and quadratic term (third model) fitted on a phase 
interval of 0.3 around the eclipse. The fourth model is as 




N (cycles) 



Fig. 1. Eclipse time delays with respect to a constant or- 
bital period model plotted vs. the orbital cycle for all the 
available eclipse time measures spanning the period from 
1977 to 2008 together with the best- fit parabola (top panel), 
and residuals in units of a with respect to the best-fit 
parabola (bottom panel). Different sy mbols indicat e differ - 
ent datasets: black squares are fr om iHellier et al.l (tT990l ). 
red stars are from IParmar et al.l (j2000l) . cvan squares are 
from RXTE data, the magenta dot is from XMM-Newton 
data, and the green triangles are from Chandra data. 

the third model plus a cubic term fitted on a phase inter- 
val of 0.4 around the eclipse. The fifth model consists of a 
Gaussian and a constant plus a sinusoid of period fixed to 
the orbital period fitted on the whole 0-1 phase interval. 
The models from the sixth one to the tenth one are as the 
fifth model, plus from 2 to 6 sinusoids with periods fixed 
to 1/2 up to 1/6 of the orbital period, respectively. The 
addition of higher harmonic components was required to 
better describe the overall orbital light curve shape, which 
differs from a pure sinusoid. We restricted our fitting to the 
first six harmonics because the addition of higher harmonic 
components was not statistically significant based on an 
F-test. 

Thus we obtained 10 eclipse arrival times (each corre- 
sponding to one of the models described above) for each 
orbital light curve. The final eclipse arrival time for each 
orbital light curve was chosen to be the average of these 
10 values, and the associated uncertainty was chosen to be 
half of the maximum range spanned by these values (1 a 
error included). The uncertainty derived in this way fully 
takes into account significant discrepancies among the dif- 
ferent eclipse arrival times found with a particular model 
to describe the eclipse and the orbital modulation. The ob- 
tained values of the eclipse epochs for each of the 10 orbital 
light curves and the relative uncertainties arc reported in 
Table m 

We then computed the eclipse time delays by sub- 
tracting from our measures the eclipse arrival times pre- 
dicted by a constant orbital period model adopting the or- 
bital period, Pnrhn, and the reference time, Tq, given by 
IParmar et al.l (|2000f ). These time delays were plotted ver- 
sus the orbital cycle number N. The integer N is the exact 
number of orbital cycles elapsed since Tq; i.e., N is the 
closest integer to (T^ — Tg)/Porbo under the assumption 
that \T^ - (Tq -I- iVPorbo)! << ^'orbo that we have veri- 



L. Burderi et al.: A Stable Orbital Period Derivative over the Last 30 Years 
Table 1. New X-ray eclipse times for 2 A 1822-371. 



3 



Ted (MJD) Error Cycle Satellite (ObsID) Tstart - Tstop (MJD) 



50352.10094 


0.00047 


20410 


RXTE (P10115) 


50352.34 - 


50353.44 


50992.0857 


0.0023 


23167 


RXTE (P30060) 


50992.81 - 


51019.65 


51975.01132 


0.00031 


27402 


RXTE (P50048) 


51975.44 - 


52101.48 


52432.03655 


0.00030 


29371 


RXTE (P70036/37) 


52432.37 - 


52435.77 


52487.97497 


0.00038 


29612 


RXTE (P70037) 


52488.36 - 


52504.01 


52519.07766 


0.00085 


29746 


RXTE (P70037) 


52519.15 - 


52547.87 


52882.09667 


0.00037 


31310 


RXTE (P70037) 


52882.02 - 


52885.21 


51975.06934 


0.00056 


27402 


XMM-Newton (230101) 


51975.55 - 


51976.15 


51779.6317 


0.0019 


26560 


Chandra (671) 


51779.70 - 


51780.16 


54607.19592 


0.00056 


38742 


Chandra (9076/9858) 


54606.96 - 


54610.52 



Note: Uncertainties are calculated as described in the text. Tstart - Tstop indicates the time interval over which a folding of the 

orbital light curve was performed to derive the time of eclipse. 



Table 2. Best-fit orbital solution for 2A 1822-371 derived 
from the analysis of the eclipse arrival times from 1977 to 
2008. 



X [(n- l/3)/(n-H5/3-0.2m2,o.i"ir^)] 



(2) 



Parameter 


Units 


P2000 


This work 


To 


MJD 


45614.80964(15) 


45614.80948(14) 


-Porb 


s 


20054.1990(43) 


20054.2056(22) 


Porb 


10-^0 s 


/s 1.78(20) 


1.499(71) 


xVd.o.f. 




21.4/16 


38.69/25 


Note: Errors 


are at 1 


a c.l. on the last 2 dU 


nts. The value of 



^orb is referred to Tg. The best-fit orbital parameters reported 
in this work are com pared with the ephemeris given by P2000 
l|Parmar et al.ll2000i ). 



fied a posteriori. These results are shown in Fig. [T] together 
with all delays computed from p reviously availabl e eclipse 
times, namely those given by He llier et al] (|1990D and by 
[Parmar et al.l (|2000) . respectively. 

These points show a clear parabolic trend that we fitted 
to the equation 

6T^ - ST^ + JPorb X iV + (l/2)Porb/'orb X , (1) 

where the correction to the adopted value of the eclipse 
time, 6Tq, and to the adopted value of the orbital period, 
^-PorbO; and the orbital period derivative, Porb, are the fit 
parameters. We get a very good fit with a ^ jd.o.j. — 
38.69/25 = 1.5. In agreement with previous results, we find 
a highly significant derivative of the orbital period, which 
indicates that the orbital period in this system is increasing 
at a rate of P^b = (1-499 ± 0.071) x 10"!° s/s. The best-fit 
values for the orbital parameters, calculated with the cor- 
rections we found from the fit of the parabolic trend of the 
eclipse epochs with Eq. [1] are shown in Table [5] Note that 
a similar orbital period derivative was recently found w ith 
new measures of optical eclipses by iBavless et al.l (|2009l) . 



3. Orbital evolution of 2A 1822-371 

Apart from mass transfer between the companion and the 
neutron star, the orbital evolution of this binary system is 
expected to be driven by the emission of gravitational waves 
and by magnetic braking. Under the further assumption of 
conservative mass transfer, orbital evolution calculations 
show that the orbital period derivative should be 



(see IDi Salvo et all (I2008I1 : IVerbuntI (|1993D : see also 
i Rappaport et al. (|1987^ )■ where mi and m are the mass 
of the primary. Mi, and the total mass, M\ -I- M2, in units 
of Mq respectively, m2, 0.1 is the mass of the secondary in 
units of 0.1 Mq, is the orbital period in units of 5 h 
(that is appropriate for 2A 1822-371 because Porb = 5.57 
h), n is the index of the mass-radius relation of the sec- 
ondary P2 oc M2 , and where the term Tmb ~ 20 takes into 

account the effe ct of the magnetic brakin g. 

In line withlVerbunt fc ZwaanI (Il98lll . Fv^W^ (|1993D . 
and iKind (|l988h (see iTaurid (|200lD for a review), we can 
parametrise this term as 



MB 



19.3 (//fco.277)"'™2^o.l"^l'^'-Pf 



2 



(3) 



where / is a dimensionless p arameter of order of unity: pre- 
ferred values are / = 0.73 (|Skumanichlll972D or / = 1.78 
(jSmithl [l979ll . and fco.277 is the radius of gyration of the 
star fc in units of 0.277, which is the appropriat e value 
for a 1 Mq ZAMS star (IClaret fc Gimened Il989ll. Note 
that the expression for the MB term given in IVerbun^ 
(|1993l ) is recovered from the abo ve adopting / = 1 and 
fc for a 1 Mq ZAMS star. Actually. iTaurij (|200ll) discussed 
three different expression for the angular momentum losses 
due to MB, Jmr, nam ely that proposed by Sku r nanic h 
()Verbunt fc Zwaanl ll981'). that proposed bv' StepienI (|1995D . 
and that proposed byfRappaport et al. (198^1), respectively. 

0.1| JskumanichI with |jRappaport| SOmC- 



However | Jg 



tcpicn 



< 



Po, 



-3.0 X 10"" 



mx TO2, 0.1 "I 



l/3p-,^/^x[1.0-fTMB] 



what between them. Because we found that to describe 
2A 1822-371 a quite large Jmb is required, we decided to 
adopt I JskumanichI which rcsulted in the term Tmb ^ 20 
adopted above. 

The orbital period derivative we measured cannot be 
explained by a conservative scenario however. A positive 
orbital period derivative certainly indicates a mass-radius 
index n < 1/3; this is indeed a quite general result, which 
does not depend on the details of the angular momentum 
losses (see also Eq. 01 below). However, the orbital period 
derivative we measured, Porb = 1-50(7) x 10~^° s/s, is about 
three orders of magnitude larger than what is expected even 
including the (strongest) MB term! This discrepancy is em- 
barrassingly large suggesting that the conservative evolu- 
tionary scenario cannot be applied in th i s case . A similar 
conclusion was reached by IBavless et al.l (|2009D , who give 
an improved ephemeris for this source based on new optical 
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eclipse measures; these authors also note that an extremely 
high mass accretion rate onto the neutron star, correspond- 
ing to about four times the Eddington limit, would be re- 
quired to explain the observed large orbital period deriva- 
tive, and conclude that much of the transferred mass must 
be lost from the system. Below we show how the orbital 
period derivative we measured can be used to constrain the 
mass transfer in the system, and how this strongly indicates 
that a large fraction of the mass which the companion tries 
to transfer to the neutron star is lost by the system. 

The mass-loss rate from the secondary can be easily 
calculated as a function of the orbital period of the system 
and the measured orbital period derivative combining the 
third Kepler law, which must be always satisfied by the 
orbital parameters of the system, with the condition that 
in this persistent system the neutron star is accreting mass 
through Roche Lobe overflow. This means that the radius of 
the secondary follows the evolution of the secondary Roche 
Lobe radius: Rl2/Rl2 — R2IR-2-, where for the secondary 
we adopted a mass-radius relation i?2 oc M2 and for the ra- 
diusof the secondary Roche Lobe we adopted the Paczvhskii 
(|1971D approximation Rl2 = 2/3'*/3[g/(l + q)Y/^a, where 
a is the orbital separation, which is valid for small mass 
ratios, q — M2/M1 < 0.8. From these conditions it is pos- 
sible to derive a relation between the mass-loss rate from 
the secondary and the orbital period derivative 



TO_8 = 3.5 X (3n — 1) ^m2, 0.1 



P-^ 



Pi 



51i 



(4) 



where rh-g is the secondary mass derivative (negative since 
the secondary star looses mass) in units of 1Q~^Mq yr~^, 
and P_io is the orbital period derivative in units of 10~^°. 
We stress that in Eq. ^ an expression for the angular 
momentum losses mechanism that drives the evolution of 
the system (e.g. MB or GR) does not explicitly appear. 
This is quite relevant because at present there is no general 
consensus on the absolute strength of the MB term nor on 
its functional dependence on the other orbital parameters. 
Indeed the effects of the driving mechanism are implicitly 
considered through the orbital period derivative, which is 
a measured quantity in our case. 

Equation (U can be inverted to derive the mass-transfer 
timescale t^^ — M2/(— M2) 



•m 



X (l-3?i)(P5h/P-io) X lOSr. 



On this short time-scale the response of the secondary star 
must be adiabatic. For m2 ~ 0.3 (see below) the envelope is 
convective and the appropriate index is = — 1/3, in agree- 
ment with the condition n < 1/3 discussed above. With this 
value we find Tj^y ~ 4 x 10^ yr for 2 A 1822-371. We note 
that the Eddington limit (in units of IO^^Mq yr^^) for ac- 
cretion onto a neutron star is ms -s — 1-54 i?g(mi), where 
i?e(mi) is the neutron star radius in units of 10^ cm, which 
slightly depends on the neutron star mass once an equation 
of state (EoS) for the ultradense matter is adopted. 

Thus, adopting m2, 0.1 > 3.3 (j Jonker et al.l l2003l ) and 
n = —1/3 in Eq. (U), we have to conclude that the sec- 
ondary mass loss rate in 2A 1822-371 is super-Eddington. 
We are therefore forced to conclude that the evolution of 
the system is highly non-conservative. 



In order to search for a possible evolutionary scenario 
for 2A 1822-371 we make the assumption that the neu- 
tron star is accreting at the maximum possible rate, i.e. the 
Eddington limit. It has to be noted that the Eddington limit 
strictly holds for a spherical geometry, and may not be a 
constraint for highly magnetised neutron stars for which the 
accreting matter is channeled onto the magnetic polar caps 
and the geometry of the matter distributi on over the Alfven 
surfac e may not be symmetric (see e.g. iBasko fc Sunvaevl 
(|l976f l). However, our assumption is justified because the 
luminosity function for highly magnetized neutron stars 
(usually found in High Mass X-ray Binaries) does not dis- 
agree with this assumption (see e.g. iGrimm et al] (p0 02f)). 
In particular, no highly magnetized neutron star is known 
to accrete at a rate much higher than the Eddington limit, 
and the most luminous high mass X-ray binaries containing 
a neutron star in our Galaxy reach luminosities of the order 
of the Eddington limit. Moreover, the extrapolated X-ray 
luminosity of 2A 1822-371 does not indicate an extremely 
high X-ray luminosity. Hence we do not have any evidence 
that the limiting mass accretion rate in this source is very 
different from the Eddington limit. 

This results in the following condition: 



(3 X m_8 — 1-54 Re{mi), 



(6) 



where {3 is the fraction of the mass lost by the secondary, 
which is accreted by the neutron star, namely Mi = —/3M2, 
where Mi and M2 are the mass derivatives of the pri- 
mary and the secondary, respectively. We consider two 
EoS, namely the moderately soft FPS and the stiffer L 



(jCook et al. IIQQJ) . which give the relation RQ(rni) for the 
neutron star radius. With this and considering the mass 
function of the system der ived from the timing analy - 
sis of the neutron star spin (|Jonker fc van der Klisl[200lh . 
Eqs. (g]) and (O can be solved to derive M2 and (3 for any 
value of TOi between mi min = 0.97 and mi max, which 
depends on the particular EoS adopted, which is 1.8 and 
2.5 for EoS FPS and L, respectively. In Fig.[2]M2, in units 
of the Eddington mass transfer rate, and /? are plotted for 
the appropriate range of neutron star masses for the FPS 
case (for the L EoS the value of /? is 30% higher, while M2 

in units of the Eddington limi t is 30% lower) . 

We now consider Eq. (3) of lDi Salvo et all ()2008[ ). which 
expresses the conservation of the angular momentum of the 
system 



(5) -Porb 



= 3 



J 

Jorh 



Rorh 

where, in this case 
j 



M2 
M, 



(7) 



— = -5.5 X 10-iS[1.0 + Tmb] mi m2, 0.1 m-^^^'P'^^/^iS) 

represents all the possible losses of angular momentum from 
the system caused by MB and GR, where Tmb is given by 
Eq.Iland 



j{(3, q,a)^l-(3q~{l~P){a + g/3)/(l + q) 



(9) 



takes into account the effects of angular momentum losses 
because of mass loss from the system, a is the specific angu- 
lar momentum of the mass leaving the system, 4j , in units 
of the specific angular momentum of the secondary, that is: 
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M 

Mem 



a 0.86 



0.82 




1.1 1.2 1.3 1.4 1.5 1.6 1.7 1. 

Ml /Mo 

Fig. 2. Secondary mass loss rate in units of the Eddington 
limit for the FPS EoS (top panel, solid line), fraction /3 of 
the mass lost by the secondary star which is accreted onto 
the neutron star (top panel, dashed line) , and specific angu- 
lar momentum of the mass leaving the system, a, in units 
of the specific angular momentum of the secondary com- 
puted for two different pres cription of the ma gnetic braking 
strength, na mely / = 0.73 (|Skumanichlll972[ ) and / = 1.78 
(jSmith 1979) (bottom panel). The vertical line at 0.97 Mq 
re presents the lower lim it of the neutron star mass derived 
bv lJonker et~al] (|2003( ). 

a = lej/{^o,hrl) = le,Porb{Mi + M2Y / {2110^ Ml) , where 
r2 is the distance of the secondary star from the center of 
mass of the system, and a is the orbital separation. 

Adopti ng the two values of / discussed a.bove, namel y 
/ = 0.73 (|Skumanichl [l972D or / = 1.78 (|Smithl I1979D . 
and ko,277 = 1, Eqs. ©, ©, (P, and © can be solved to 
derive a as a function of mi. In Fig. [2] a is plotted for the 
appropriate range of neutron star masses for the FPS case 
(for the L EoS the values of a are 7% higher). 

The values of a we obtain are in between the spe- 
cific angular momentum at the inner Lagrangian point, 
= [1- 0.462(1 + g)2/3gi/3]2 ^0.4 for 2A 1822-371, and 
the specific angular momentum of the secondary, a2 = 1-0, 
and actually quite close to a2 ■ This is expected if the mass 
lost by the secondary star is blown away because of the ra- 
diation pressure exerted by the Eddington luminosity gen- 
erated by the accretion onto the neutron star. 

For both the adopted EoS and all the possible values 
of the neutron star mass, the values of /3 are in the range 
0.13 — 0.29, which means that the mass transfer in 2A 1822- 
371 is not conservative, at least, at 70% level, which, as 
we already noted, is true independently of any assumption 
on the particular angular momentum losses . Interestingly, 
this is the key that opens the possibility of constructing 
a consistent secular evolution for this system. Indeed the 
contact condition, Rl2/Rl2 = R2/R2, can be solved to 
derive a theoretical prediction for the mass-loss rate once 
a prescription is given for the possible losses of angular 
momentum from the system caused by MB and GR as in 
Eq. (©: 



-3.5 X 10-''[1.0 + Tmb] 



TOi m2 I rn' 



F{n,giP,q,a)), 
where 

F(n,g(/3,g,a)) = [n - 1/3 + 2gip, q, a)]-\ 



(10) 
(11) 



The function F(n, g(/3, 5, a)) is very sensitive to the sce- 
nario adopted: for a conservative scenario (/? = 1) 
F{n, g{/3,q,a)) ^ 1.5 while for (3 in the range 0.13 - 0.22 
(which is appropriate for the FPS EoS and all the possi- 
ble values of the neutron star mass) F{n,g{P,q,a)) ~ 40. 
This means that the term acting to shrink the secondary 
Roche Lobe - because of the extra angular momentum 
losses caused by the mass expelled from the system - de- 
termines an amplification of the mass-loss rate through the 
function F(n, 5(/3, g, a)) with respect to a fully conservative 
case. This amplified mass-loss rate in turn determines the 
high value of the orbital period derivative observed in this 
system. Consequential angular momentum losses (CAML), 
i.e. angular momentum losses that are themselves the result 
of mass transfer, have been proposed in the c ontext of cat- 
aclysm ic variables (CVs) evolution (see e.g. iKing fc Kolbl 

Inserting the values determined in this paper for 
2A 1822-371 in Eq. 1^ we find r^-^ = 2.1 x 10^ yr. This 
means that the system as it is observed now will proba- 
bly end on this timescale, possibly with the tidal disrup- 
tion of the companion star. Indeed, iKing fc Kolbl ([l99l) 
argued that mass transfer could be unstable when CAML 
are present. This time-scale is extremely short, which in- 
dicates that it is possible that some short orbital period 
LMXBs can last much shorter than what was previously 
thought. This evolutionary phase, characterised by a super- 
Eddington mass transfer rate, may be a common phase in 
the evolution of LMXBs, albeit short-living. Because this 
phase should not last more than a few million years, there 
may be very few observed systems in this phase (e.g. the 
so-called Z-sources, which are persistently bright LMXBs). 
This could have profound implications for the estimate of 
the actual number of LMXBs produced in the Galaxy as 
inferred from the observed ones, and also for the predicted 
number of millisecond binary pulsar. We note that this 
would help to bring the number of LMXBs in line with 
the estimated number of millisecond binary pulsars. But a 
detailed analysis of this delicate and long-standing prob- 
lem needs a dedicated study of this almost unstable phase 
of the orbital evolution, which is beyond the scope of this 
paper and will be discussed in a forthcoming paper. 
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